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Abstract 

Propagation of a particle accelerated by an external field through a scattering 
medium is studied within the generalized Lorentz model allowing inelastic 
collisions. Energy losses at collisions are proportional to (1 — a 2 ), where 
< a < 1 is the restitution coefficient. For a = 1 (elastic collisions) there 
is no stationary state. It is proved in one dimension that when a < 1 the 
stationary state exists . The corresponding velocity distribution changes from 
a highly asymmetric half-gaussian (a = 0) to an asymptotically symmetric 
distribution ~ exp[— (1 — a)v i /2], for a — ► 1. The identical scaling behavior in 
the limit of weak inelasticity is derived in three dimensions by a self-consistent 
perturbation analysis, in accordance with the behavior of rigorously evaluated 
moments. The dependence on the external field scales out in any dimension, 
predicting in particular the stationary current to be proportional to the square 
root of the external acceleration. 

PACS numbers: 51.10.+y 
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I. INTRODUCTION 

The object of the present paper is to study the propagation of a particle through a 
medium composed of immobile spherically symmetric scatterers (of infinite mass), randomly 
distributed in space with some number density n. Between collisions the particle moves with 
acceleration a, acted upon by a constant and uniform external field. At binary collisions its 
velocity v is instantaneously transformed according to the law 

v — ► v' = v — (1 + a)(v • a) a (1) 

where < a < 1 is the so called restitution coefficient, and a denotes the unit vector along 
the line passing through the centers of colliding particles at the moment of impact. Whereas 
the projection of velocity v on the direction tangent to the surface of the scatterer at the 
point of contact is not changed, the component (v • a) along the normal gets multiplied by 
(—a). In the extreme case of a = 0, the postcollisional velocity reduces to the tangential 
component. 

When a — 1, collisions ([[]) are perfectly elastic, and we recover the situation of the 

classical Lorentz model of electric conductivity in metals [Jj. It can be proved that there is 

no stationary state in this case 11,0. This is because the moving particle does not loose 

energy at collisions, and gets asymptotically heated by the field beyond any bounds. In 

particular, its mean kinetic energy diverges as ~ t 2 ^ 3 when the time t — > oo. 

The aim of this work is to show, that in the case of inelastic collisions (H) the stationary 
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state becomes possible in the whole range of < a < 1. Once a < 1, there occurs 
dissipation of the kinetic energy E = mv 2 /2 at encounters, as in accordance with (|1|) E 
suffers the transformation 

E -> E' = E - m(l - a 2 )(v ■ a) 2 /2 (2) 

The inelastic dissipation (0) turns out to suffice to balance the energy flow from the external 
field. 

The basis for the subsequent analysis is the linear Boltzmann equation satisfied by the 
probability density /(v; t) for finding the propagating particle with velocity v at time t. In 
writing this equation one must take into account that the jacobian of transformation ([!]) 
equals a (dissipative collisions are contracting the volume in the velocity space), and that 
the inverse transformation is obtained by replacing a in ([!]) by a" 1 . The kinetic equation 
reads 

^ + a • /( V ; t) = M J da(v ■ a)6{v • *) {^V[v - (1 + a" x (v ■ *)*); t] - /(v; t)} 

(3) 

Here A = (nnR 2 )^ 1 denotes the mean free path (R is the sum of the particle and the scatterer 

radii), v is the unit velocity vector, and the integration with respect to da spreads over a 

unit sphere (solid angle). In the gain term (due to collisions) the factor a' 2 compensates 

for the contraction of the velocity space and gives the proper value to the precollisional 

velocity of approach (v • a) /a. Owing to this factor the velocity integral of the collision 
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term vanishes, which permits to deduce from ([3]) the continuity equation. 
A convenient equivalent form of equation (|3]) reads 

f 9 d \ „, . Ivlf 9 r dn „ rl + a. . „ 1 — a 1 . , .1 . . 

[m + a ' ef) /(v; () = T K / -J [— |v|n - — v; ( J - /(v; () | (4) 

It can be obtained from (Q) by changing the angular integration variables from da = 
sinipdipdrf) (where cosip = (v • a)) to dn = sinxdxd<p, with x = 2^. When a = 1, we 
recover the kinetic equation derived originally by Lorentz Q, 

Our work is closely related to the study of spontaneous percolation in three dimensions 
of Wilkinson and Edwards @. The dynamics of particles falling under gravity through a 
fixed scattering medium was analyzed therein on the basis of a Boltzmann-like equation 
with a phenomenological scattering function. The perturbative methods developed in 
permitted to analyze the low inelasticity limit leading to the same qualitative predictions 
as those derived here. The novelty of our contribution consists mainly in proving by an ex- 
plicit construction the existence of a stationary state in one dimension, and in developing an 
original self-consistent perturbative expansion in three dimensions, inspired by the scaling 
structure appearing in one dimension. It is to be stressed that the rigorous one-dimensional 
results firmly support the predictions of the perturbative analysis in three dimensions, as 
the properties of the system in both cases are qualitatively the same. We could also cal- 
culate rigorously some of the moments of the stationary distribution in three dimensions. 
Their asymptotic behavior in the elastic limit coincides with that predicted by the proposed 
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perturbation scheme. In Section II the dimensional analysis of the stationary solution to the 
kinetic equation (^) is performed. Then, in Section III the one- dimensional version of the 
model is considered. We prove therein the existence of a stationary solution and we analyze 
its dependence on parameter a. Section IV is devoted to the extension of our results to 
three dimensions. This is obtained in the low inelasticity limit a — > 1 by a self-consistent 
construction of a stationary state with appropriate scaling behavior. The paper ends with 
concluding comments. 

II. DIMENSIONAL ANALYSIS 

It is as usual instructive to formulate the theory in terms of dimensionless variables. We 
thus put 

v = VaXu (5) 

where a = |a|. Denoting by F(u) the stationary dimensionless distribution we find from (|J) 
the equation 

9 s , ( 9 fdri rl + a. _ 1 — a 1 „. .1 
a • —F(u) = u I a~ 2 / —F u n u - F(u) } 6 

The amplitude a of the external field disappeared from equation (H). Only the unit vector 
a shows there, indicating the direction of acceleration. Hence, the field dependence of the 
moments of distribution F can be obtained by the velocity scaling (|5|). In particular 

< v >= VaX < u >~ y/a, < v 2 >~ a (7) 
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The fact that the stationary particle current is proportional to the square root of the field 
is of no surprise. The immobile scatterers represent a zero-temperature medium, and the 
only energy scale comes from the acceleration a. The relations (|?|) could be thus predicted 
on purely dimensional grounds. At fixed restitution coefficient a, one cannot expect the 
linear response even for very weak fields. 

The above conclusions apply in any dimension. In the study of the stationary state we 
shall use in the coming sections the dimensionless form @ of the kinetic equation. 

III. ONE-DIMENSIONAL STATIONARY VELOCITY DISTRIBUTION 

In one dimension equation (Q) takes the form 

(8) 
Putting 

F(u) = G(u\u\/2) (9) 

we find 

G'{s) + G{s) =a~ 2 G{-a- 1 s), (10) 
where G' denotes the derivative of G. Applying to (PUD the Fourier transformation one finds 

oo 

G(k) = (1 + iky^-afk) = G(0) nt 1 + K-^Yk]- 1 (11) 

r=0 

where 
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n n i 



F(u) 



u 



a 



2 F(-a 



u 



F(u) 



G(k) = J 



dse- iks G(s) 



For < a < 1 the infinite product converges which proves the existence of the stationary 
state. This is an important conclusion indeed, showing the fundamental difference with re- 
spect to the Lorentz model with elastic collisions. Clearly, an arbitrary degree of inelasticity 
suffices to balance the energy flow from the external field. 

A particularly simple result is found for perfectly inelastic collisions. Indeed, when a = 
equation (]TTf) reduces to 



Inverting the Fourier transformation and using (|9]) we find 



In one dimension each collision with a = completely stops the particle, dissipating the 
whole energy absorbed from the field. So, in a stationary flow the velocities oriented against 
the field are not possible (the 9 factor in (0)). It is quite remarkable that the half space of 
possible velocities gets a gaussian weight. 




(13) 



Let us turn now to equation fllPl) . It is equivalent to the system of two coupled equations 



G' + (s) + G4s) 



a- 2 G_(s/a 2 ) 



(14) 



G'_(s) + G+{s) = a- 2 G+{s/a 2 ) 
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where G + (s) = [G(s) + G(—s)]/2 and G_(s) = [G(s) — G(—s)]/2 are the symmetric and 
antisymmetric parts of G, respectively. The second relation in (|TJ]) can be solved for G_(s) 
yielding 



GL(s) = y cfeG+(z) (15) 
We then get from (0) a closed equation for the symmetric part 

— 2 —4 

/•sq /'-sa 

G' + (s) = - dzG + {z)-a- 2 ^ dzG+(z) (16) 

J set 2 

In order to study the low inelasticity limit we define a small parameter e = 1 — a. When 
e< 1, equation (1TB) takes the simple asymptotic form G'As) = — 4esG + (s), and we find 



G+(s) = Ce 1/4 exp(-2es 2 ) (17) 



Using then the relations (15) and (^]) we can determine the asymptotic form of the velocity 



distribution for weakly inelastic scattering. It reads 

Ce 1/4 

F(u) = ^^{l + eu\u\)exp(-eu 4 /2), e<l (18) 
v2 

The value of constant C is independent of e and follows form the normalization condition 

C- 1 = J dwexp{-2w 4 ) 

The dominant term in (|18|) is symmetric and depends on the scaled variable w = eu A . It 
implies the divergence ~ e -1 / 2 of the kinetic energy < u 2 > for e — ► 0. This divergence 
reflects the disappearance of the stationary state in the case of elastic collisions. In fact, 



the kinetic equation (H) permits to calculate directly some of the moments of distribution 
F, and to predict divergences in the elastic limit. An interesting example is provided by the 
recurrence relation 

2j — 1 

<u j \u\ j >= ; rr— — < v?~ x \u\° >, with < u\u\ >= (1 + a)' 1 (19) 

11 1 - (-lya 2 !- 1 11 v ' 

which follows directly from when mutiplied by u^\u\^~ l and integrated over the velocity 
space. For j = 2, we get from (TBI) 
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< u >= 3[(l + a)(l - a 3 )] 

implying the divergence ~ e -1 of the fourth moment. It should be stressed that the asymp- 
totic state (0) is consistent with the exact relations ([19]). This result can also be obtained 
directly from the explicit solution by studying the asymptotics of the infinite product in 
(|TT|) as a. — > 1. The high kinetic energy for e — > implies a large collision frequency, which 
makes the asymptotic distribution (|T8|) almost symmetric. In this situation the stationary 
current is very weak. Indeed, we find 

<u>=Ct 1/A /V2 (20) 

Before closing let us note that the above results can be generalized to the case where the 
collision frequency is proportional to |w| 7 , < 7 < 1, rather than to \u\. The corresponding 
generalization of the kinetic equation (H) reads 

F'(u) = \up la-^Fi-a^u) - F{u)} (21) 



and the substitution 

F(u) = G[u|w|V(7 + 1)] 

transforms (pT]) into 

G'(s) = - G(s) (22) 

Performing then the same analysis as for 7 = 1, we find the generalization of the asymptotic 
formula (|18|) 

F(u) = Ce 1 / 2 ^ +1 )( 7 + l)-V7+i(i + eu|uDexp[-en 2 ^ +1 V(7 + 1)], e < 1 (23) 

Again the dominant term in (^) is symmetric, and the first two moments follow the 
asymptotic law 

<it>~e^ +1 >, <u 2 >~e-^ +1 

An interesting case is that of the so called Maxwell gas, where 7 = 0, so that the collision 
frequency becomes independent of the velocity. The stationary current stays then finite even 
in the limit e — > 0. However, the ~ e _1 divergence in the kinetic energy indicates the infinite 
absorption of energy from the field for e = 0. 

Our analysis for 7=1 (free motion between collisions) shows that the stationary velocity 
distribution changes from a highly asymmetric half-gaussian (f[3|) at the strongest dissipation, 
to an almost symmetric scaled distribution ~ exp(— ew 4 ) in the weak dissipation limit. The 
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experience from the studies of the classical Lorentz model showed the independence from 
the spatial dimension in the qualitative behavior. We can thus expect the appearance of the 



same scaling structure in three dimensions. 



IV. STATIONARY VELOCITY DISTRIBUTION IN THREE DIMENSIONS 



In three dimensions the stationary velocity distribution F(u) satisfies the equation (see 



!)) 



u < a 



2 I dA F 

47T 



1 + a, 



a 



2a 



u n 



2a 



-u 



F(u) 



(24) 



Multiplying ( f24|) by a function vjj(u) and integrating over the velocity space one obtains 



the dual equation of the form 



1 — a 1 + a . 



-u 



u n 



- U tfl U > > 



(25) 



Here < ... > denotes the average with respect to the probability density F. One can obtain 
from (E3) useful information by an appropriate choice of tp. In fact, some moments of F can 
be evaluated exactly in this way for any value of a. So, for instance putting ip(u) — u • a 



yields the formula 



1 + a) < |u|(u ■ a) >= 2 



(26) 



And when ip = |u| 3 , we find 



3 < lulfu ■ a) >= 



2(1 - a 5 ) 
5(1 -a 2 ) 



< u > 



(27) 
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Combining (|26"D and ([27]) permits to determine the fourth moment 

(1 - a)(3 + 6a + 4a 2 + 2a 3 ) < |u| 4 >= 30 (28) 

In accordance with the remark made at the end of section III, we find here the same di- 
vergence < |u| 4 >~ e _1 for e = (1 — a) —>■ as in one dimension. Let us finally note the 
relation 

4 < (u-a) >=< |u| 3 > (1 -a 2 ) (29) 

following from ( p5|) for i/)(u) = |u| 2 . 

In order to determine the distribution F in the limit of weak inelasticity we put a = (1— e) 
in equation (p4|), and expand the inelastic collision law up to terms linear in e. Clearly, the 
stationary state can depend on two variables only, u = |u| and \x = a • u. Expressing the 
differential operator in (|24j ) in terms of u and fi one finds eventually the following asymptotic 
equation 

(4 + ! + ") F4u - = -"s Fs < m ) + f (iiif" 4Fs Mi - -Jy /-, *"^(«.»)]) 

(30) 

By analogy with the approach developed in the one- dimensional case, we introduced here 
the spherically symmetric projection of F 

Fs{u) = J d £F(u) 

and the deviation form spherical symmetry F& = F — F$. 
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An important relation between F$ and Fa is obtained by integration of fl3~0|) over the 
whole range of the angular variable —1 < fi < 1. One finds 

J d/j, hFa(u, n) = eu 2 F s (u) (31) 

In order to proceed, we shall assume that for e < 1 the state F depends on velocity 
via a scaled variable e x u. The exact relations fl2T)|)-(P§D imply then the value x = 1/4. 
Indeed, we know from ( PB[ ) that the fourth moment diverges like e _1 , whereas the assumed 
scaling predicts the ~ e~ 4x behavior. So, taking also into account relation (|5ID , we write the 
perturbative expansions of Fs and Fa as 

FJH = F° s (w) + VeF^w) + eF 2 s (w) + ... (32) 
F A (w, n) = VeF A (w, fj) + eFl(w, fj) + 

where w = e 1//4 w. 

The expansion of the non-spherical part Fa does not contain the zero order term in ac- 
cordance with the relation ([H|) between Fs and Fa- Inserting expansions (j3^) into equations 
(p0|) and ( |3~TD to lowest order in e we find 

wF\(w,fi) = -^Fl(w), J^dtifiF\(w,tj) = w 2 Fl(w) (33) 

From (^3j) there follows a closed equation for Fg 

3w 3 Fl(w) = -2-^F° s (w) 
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The lowest order terms in expansion (^) are then readily determined. The asymptotic 
formula for the stationary distribution reads 

F(w, n) = F$(w) + y/eF\(w, ll) = Ce 3/4 (1 + 3y/eLiw 2 /2)exp(-3w 4 /8), (34) 

where the constant C assures the normalization. Inserting into fl3~4|) w = e 1//4 w and ll = a- u, 



we obtain the same structure of the state F as in one dimension (compare with (0)). The 
same e-dependence of the moments of F thus also holds. 

This perturbative analysis can be continued. We checked that at the next step the 
following relations are found 

F s( w ) = 0, 

F 2 A (w^) = PM[l + w^-]Fl(w), 

where ^(/-O = (1 — 3/i 2 )/2 is the second Legendre polynomial. This confirms self-consistency 
of expansion (]32f) in powers of y/e. 

V. CONCLUDING COMMENTS 

We proved the existence of a stationary state in the case of a one-dimensional propagation 
of a uniformly accelerated particle suffering inelastic collisions with randomly distributed 
scatterers. The velocity distribution, strongly asymmetric at high dissipation, showed the 
scaling behavior in the elastic limit e = 1 — a — > 0, recovering asymptotically the spherical 
symmetry F ~ exp(— eu 4 ). 
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The same scaling properties could be derived in three dimensions by a suitable pertur- 
bation method. Of course, it would be desirable to provide a rigorous existence proof also 
in three dimensions. 

In the one-dimensional case we could determine the form of the velocity distribution for 
perfectly inelastic collisions (see (|13|)). In three dimensions, passing in the dual equation 
(|5|) to limit a — > 0, one can deduce in a straightforward way the form of the equation 
satisfied by the state F at a = 

vra ■ 7p^( u ) = J dw 5(w ■ u)F(w + u) - tt|u|F(u) (35) 

It should be realized that in contradistiction to the one-dimensional case the particle is 
not stopped at encounters with a = 0, but its postcollisional velocity reduces then to the 



component tangent to the surface of the scatterer. It is thus clear that the solution to (p5|) 
can give non-zero probability weight only when a ■ u > 0. Solving equation (|35D is another 
open problem related to our study. 

Finally, an interesting question would be to explore the dynamics of approach to the 
stationary distribution. This problem has been solved only in one dimension for a = ||. 
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